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ABSTRACT 


This report presents a treatment of the various factors involved in Ildar 
data acquisition and analysis. This treatment highlights sources of fundamen- 
tal, systematic, modeling, and calibration errors that may affect the accurate 
interpretation and calibration of lidar aerosol backscatter data. The discus- 
sion primarily pertains to ground-based, pulsed-CC^ lidars that probe the 
troposphere and are calibrated using large, hard calibration targets. However, 
a large part of the analysis is relevant to othet types of lidar systems such 
as lidars operating at other wavelengths; continuous-wave (cw) lidars; lidars 
operating in other regions of the atmosphere; lidars measuring nonaerosol 
elastic or Inelastic backscatter; airborne or Earth-orbiting lidar platforms; 
and lidars employing combinations of the above characteristics. 
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SECTION 1 


INTRODUCTION 

Measurements of the aerosol volume backscatter coefficients at 
CO 2 laser wavelengths not only provide an improved understanding of large 
aerosol particles in the atmosphere but also are very important for assessing 
the requirements of Earth-orbiting CC^ lidars for measurements of atmospheric 
wi .ds and other parameters. In particular, remote measurements of winds on a 
global scale are essential for a better understanding of processes taking 
place within the Earth's troposphere. In order for us to understand climatic 
and weather-related processes, measurements must be made of the wind vectors 
as a function of latitude, longitude, altitude, time, and weather conditions. 
Studies indicate such tropospheric-wind field data would greatly enhance the 
accuracy of the National Weather Service's 12 hour-5 day forecasts. There is 
general agreement that the infrared coherent Doppler lidar technique, of all 
the candidate passive and active atmospheric measurement techniques that may 
be implemented from an orbiting platform, offers the greatest potential for 
global tropospheric-wind measurements. Since the return signal from a given 
range is directly proportional to the aerosol backscatter coefficient in that 
region, uncertaintier in the knowledge of the backscatter coefficient imply 
corresponding uncertainties in the required laser-transmitter pulse energy and 
receiver-telescope diameter. Therefore, measurements of the aerosol backscat- 
ter coefficient at CO 2 laser wavelengths as a function of latitude, longi- 
tude, altitude, time, and weather conditions are critical. 

A number of CC^ lidar systems are in use or are proposed that 
will contribute to enlarging the aerosol backscatter coefficient data base. 
These include fixed-location, ground-mobile, and airborne systems; and both 
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pulsed and continuous-wave (cw) focused. At JPL, we have been operating a 
fixed, pulsed-CC^ lidar over the past two years in order to measure the 
tropospheric aerosol backscatter coefficient in the 9 to 11 urn region. 

Since quantitative values of the backscatrer coefficient are 
desired, a technique for accurately calibrating the backscatter data has been 
developed that uses a large, hard surface as a Ildar calibration target. The 
use of the total backscatter to Rayleigh backscatter ratio technique, which is 
often used in the visible wavelength region to calibrate Ildar responsiveness 
and consequently to obtain aerosol backscatter coefficients (Russell, Swissler, 
and McCormick, 1979), does not apply at mid-infrared wavelengths. The 
assumption that backscatter and extinction are related through a functional 
relationship, such as a power-law relationship (Klett, 1981, 1983), also does 
not apply. In the mid-infrared the Rayleigh backscatter is extremely low, and 
the atmospheric extinction is usually dominated by molecular species rather 
than by the aerosol itself. Thus a hard-target calibration technique is used. 

The measurement of a single vertical profile of aerosol backscatter 
requires many individual shots of the lidar system in close temporal proximity 
both in the vertical direction through the atmosphere and at the calibration 
target. Each lidar pulse results in a digitized record of transmitted power 
vs. time and received (backscattered) power vs. time. These data records must 
be processed appropriately and then combined with the pertinent characteristics 
of the lidar system and parameters of the atmosphere and calibration target in 
order to calculate the desired aerosol-backscatter range profile correctly. 

The purpose of this report is to discuss many aspects of this 
complicated and multifaceted calibration and data interpretation process. It 
is hoped that this discussion will contribute to an accurate and consistent 


calibration of Ildar data* Although he authors' primary experience is with a 
stationary, pulsed-CQ 9 Ildar system, much of the discussion will be useful 
in calibrating backscatter data from ground-mobile and airborne lidar systems, 
from Ildar systems operating at other wavelengths, and from cw-focused lidar 
systems. 
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SECTION 2 


LIDAR BACKSCATTER SIGNALS 
2.1 ATMOSPHERIC AEROSOL BACKSCATTER SIGNAL 

Consider a pulsed infrared lidar system directed into the 
atmosphere. Assume colltnear or near-coixinear (side-by-side) geometry for 
the transmitter and receiver telescopes and assume that both telescopes 
subtend very small solid angles at the ranges of interest. Let the 
transmitted pulse start at time t = 0 and end at time t * T with a power 
profile of P^Ct) [W ] , as shown in Figure la. The received power [W] due to 
aerosol backs cat ter is then given by 


p b ( t ) 


ct/2 

/ 


P tb (t- f > MR) 


c(t~T )/2 
p 


R 

~ n 0(R) exp{-2 fa. (R')dR'}dR 
R^ 0 b 


(1) 


where A is the effective receiver area, R is the llne-of-sight range from the 
_2 

telescopes, AR = to^(R) [sr] is the solid angle subtended by the receiver, 

r) is the system's optical efficiency, 0(R) is the range-dependent telescope 

overlap function defined as the fraction of the transmitted pulse energy that 

is within the receiver's f ield-of-view (FOV), [m is the total 

extinction coefficient of the atmosphere along the optical path, 

2 "“3 "*1 — i _ j 

£(R) [mm sr = m sr ] is the aerosol volume backscatter coefficient 
defined as the fraction of incident energy scattered in the backward direction 
per unit solid angle per unit atmospheric length, and where the integration 
over R indicates that the received power at time t is due to contributions from 
a slab of atmosphere of thickness ct /2 centered at R, = c(t/2~r /4). (A list 
of symbols used in this and following equations, along with their definitions, 
can be found in Appendix A.) Other sources of received power at time t, such 
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Figure 1. Pulsed lidar temporal profile examples for 

(a) the transmitted laser pulse power, (b) the 
backscattered power from the atmospheric aerosol 
and (c) the backscattered power from a hard 
tareet 






os multiple scattering, Rayleigh scattering, resonant fluorescence, fluores- 
cence, and Raman scattering, are neglected. Monochromatic transmitted radiation 
is assumed, and thus the wavelength dependence of p(R), n, 0(R), and a b (R) 
is understood. If the receiver employs heterodyne detection, then (3(R) in 
Equation (1) refers only to the bnekseattored radiation with polarization 
parallel to the polarization of the local oscillator. Figure lb depicts the 
aerosol backscatter signal. 

The desired quantity is the backscatter coefficient range profile, 
p(R). From Equation (1) we see that the maximum spatial frequency of |3(R) 
along the optical path that can be determined is on the order of 2/(ctp). 

This causes a bias against the measurement of low values of (3(R) that might 
occur over spatial distances smaller than cx^/2. The smoothing function 
that operates on (3(R) is not rectangular but is shaped by the range-dependent 
tfrms in Equation (1). If we now assume that both (3(R) and 0(R) are slowly 
v irylng compared to the spatial distance ct^/2, then Equation (1) becomes 

ct/2 R 

P b (t) »= R(R b ) A n 0(R b ) j P tb (t- ~) R" 2 exp {-2 J & b (R # )dR # }dR (2) 

C ( t~T )/2 0 

p 

At this point we are able to solve for p(R) if all of the terms in Equation 
(2) are known accurately. Equation (2) Is still exact under the stated assump- 
tions. A further simplification of Equation (2), which is commonly used, is 
to assume that the entire Integrand except P b (t) is constant over the 
Integration range, yielding 
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( 3 ) 


P 



et/2 

» B(R b ) -§~n 0(R b ) exp {“2 f « b (R')dR'} f P tb (t- -~)dR 

R b 0 c(t-T )/2 

P 

R b T P 

» A( R b > “T“ n °(R b > exp{-2 f a b (R')dR r } § f P tb (t) dt 
R b 0 0 


This simpllf ieation can load to errors tn the calculated value of P(R) and 
will be discussed later. The remaining integral in Equation (3) is simply the 
pulse energy (J]{ 


T 

/ 

E tb’j p tb (tHt w 

0 

Often the pulse power profile of Figure la is assumed to be rectangular with 

power P,,« tais case E,^ * P n T . 

0 i;b 0 p 

2.2 HARD-TARGET BACKSCATTER SIGNAL 

Consider the same pulsed lidar system directed at a calibration 

target that is at a range R from the telescopes. Let the target be larger 

s 

than both the laser spot size at R and the receiver f ield-of-view at R , 

s s 

and let the polar angle between the optical axis and the average surface 
normal of the target be 9 . If the pulse in Figure la (P (t) in this hard- 

uS 

target case) is transmitted, the received power [W] is given by 
/ 2R \ * * 

P s (t) = P tsr“rj pt ^2“ n ° (R s ) exp{ " 2 / a s (R ' )dR ' } (5) 

s 

where a s (R) tm *] is again total extinction and is most likely different 
from « b (R) due to the different atmospheric paths employed, p* [sr *] 


7 



is the target parameter defined as the reflected power per steradian toward 
the receiver divided by the incident power (Kavaya et al., 1983), 2R f ./c is 
the round-trip transit time of light, and aerosol backscatter is added to the 
earlier list of neglected sources of received power at time t. Note that the 
received-power profile is identical to the tronsmitted-power profile, but is 
delayed by 211^/c seconds. This is depicted in Figure 1c. We will neglect 
the spread in the distance to the target caused by the target angle 9 and 
the finite beam diameter. 

2.3 SINGLE-PULSE HARD-TARGET CALIBRATION OF P(R) 

Let us assume a very simple case in which one lidar pulse is fired 
at the atmosphere and one pulse is fired at the calibration target. Although 
Equation (2) is more exact than Equation (3), seldom are all the terms in 
Equation (2) known accurately. Furthermore, often a signal proportional to 
transmitted pulse energy E ^ is recorded for normalization and, therefore, 
an expression for received power that is proportional to E^, such as 
Equation (3), is desired. Comparing Equations (3) and (5), we see that the 
integral of the target return profile given by Equation (5) must be 
calculated. This integral I [J] is given by 


I 

ft 



2R 

— 2. + T 
c p 

p (t)dt = P* -|-n 0 (R ) 

S R Z o 

s 


exp {-2 / a (R')dR') E 

n S 


( 6 ) 


(We are assuming at this point that a sufficiently fast transient digitizer is 

used to record P (t) so that an accurate integral is possible. Alternatively, 
s 
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an integrating analog-to-digital (A/D) converter with a gate width that 
captures the return pulse may be used if the signal voltage is proportional to 
P (t).) Combining Equations (3) and (6) we see that 




E t b 


ts * 2 

k ° ® 


r 2 nrR v exp {-2 f a (R')dR') 

R b ° ( V 0 b 

r 2 OCR.) R, 

s exp{-2 / a (R' )dR'} 

0 b 


(7) 


where we let the atmospheric and target pulses have energies and E tg) 
respectively. Even in this simple example in which fluctuations of the 
returned signal intensity due to speckle and atmospheric turbulence are 
Ignored, and in which an adequate slgnal-to-noise ratio is assumed from one 
lidar pulse, it is clear that the lidar pulse energies; the target parameter, 
P*; the telescope overlap function, 0(R); and the extinction profiles, 

•' v . (R) and «, (R), must be well-known. In addition, the assumptions and 
simplifications leading to Equations (3) and (6) cannot be Ignored. 


SECTION 3 


DATA ACQUISITION ERRORS 

The discussion leading to Equation (7) assumed that the transmitted- 

pulse profiles, P^t) and P tg (t); the aerosol backscatter signal, P^(t); 

and the target backscatter signal, P (t), vrere all known In their fundamental 

s 

units of watts. In reality, detectors are used to convert the optical fields 
into a voltage (or current), followed by preamplifiers, amplifiers, filters, 
etc. Finally, this voltage is recorded for later processing. We may lump all 
of the characteristics of this chain of components (or receiver) into a single 
function that operates on the optical signal power: 


V(t) = F 1 {P(t)} 


( 8 ) 


where V(t) [V] is the recorded signal available for processing and is not 
necessarily a linear function of P(t). We may factor the function Fj into 
various typical subfunctions, such as the gain G of the receiver and the 
largest time constant t of the receiver, which represents the slowest element 
in the chain. Thus Equation (8) can be written 


V ( t ) = G x F 2 {P(t)} * JjexpC-t/x) x H(t)J (9) 

where the symbol * stands for the convolution operation and H(t) is the 

Heaviside unit step function. (Convolution in the time domain corresponds to 

multiplication in the frequency domain. In Equation (9) the spectrum of 

2 - 1/2 

G x F^,{P(t) } is multiplied by [1 + (2nft) ] . ) In general, the recorded 

voltages from the transmitted pulse, the aerosol br°kscatter , and the target 
backscatter are given by 
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( 10 a) 


V t (t> - G t X F t {P t (t)} * 

~ exp(-t/x t ) X H(t)J 

V b (t) * G b x F b {P b (t)} * 

exp(-t/x b ) x H(t)J 

V (t) - G x F {P (t)} * 

7 exp(-t/x ) x H(t)l 

s s s s 

L t s s j 


( 10 b) 

( 10 c) 


Wc are assuming that the slowest element of the receiver is much slower than 

the next slowest element and that this element can be approximated by a simple 

RC low-pass filter (x = RC), and we also assume that any DC offsets that are 

added to V(t) can be determined and removed. The impulse responses of the 

low-pass filters, the exponential terms in brackets, are normalized to have an 

area of 1. Thus the area under V(t) is equal to the area under G x F{P(t)} 

(Bracewell, 1978). The centroid of V(t) occurs 'i s after the centroid of 

G x F{P(t)}. The spread of V(t) about its centroid (similar to variance) is 
2 2 

x s greater than the spread of G x F{P(t)}. Frequency components of 
G x F{P(t)} much less than 1/x are not attenuated, while frequencies much 
greater than 1 /t fall as 1/f. If x approaches 0 (as with very fast 
receiver components), then the RC-filter impulse response approaches d(t) 
and the convolutions in Equation (10) can be neglected. 

Strictly, the first four terms in the right side of Equation (7) 
must be replaced by appropriate expressions using Equation (10). Equation 
(10a) must be applied twice to solve for P , (t) and P (t), which are then 

L D L S 

substituted into Equation (4) to obtain E^ and E^. Equation (10b) must 
be solved for P, (t) and Equation (10c) must be solved for P (t), which is 

D S 

then substituted into Equation ( 6 ) to obtain I . Fortunately, several sim- 
plifications of this procedure are possible. It was stated earlier that the 
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convolution operator in Equation (10) does not affect the area under V(t), but 
does spread V ( t ) in time with respect to F{P(t)}. Since Equation (4) is used 
twice with Equation (10a) to take time integrals, the convolution in Equation 
(10a) may be neglected, provided that the upper integration limit in Equation 
(4) is extended to include all of V t (t) and provided that the function 
is linear. We will make these assumptions and thus Equation (10a) reduces to 


v tb<‘> - 




Ola) 


v ts ( t) = 


G t p t.< e > 


(lib) 


The assumption that the function F is linear is reasonable since linear 
detectors are usually used to monitor pulse power. Note in Equation (11) that 
it is also assumed that the gain, G^, is constant. Actually, G^ is a 
function of the wavelength, X, and may also depend on P and time. If only 
a small portion of the pulse cross section is Imaged onto the pulse-power 
detector, then pulse-to-pulse variations in the pulse spatial profile may 
cause variations in the fraction of pulse power incident on the detector. In 
effect, a pulse-to-pulse (or time) variation in G results. 

It is not quite as reasonable to eliminate the convolutions in 
Equations (10b) and (10c), since often a lidar receiver is not linear in 
received power. Two common nonlinear examples of the functions F^ and F^ 
are /P or AnP. The fact that P^(t) probably does not have frequency 
components of interest greater than 1 /t u allows the convolution in Equation 
(10b) to be eliminated. Even though Equation (6) is used to integrate P (t)> 
the convolution in Equation (10c) may not be ignored unless F g is linear or 
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t , Is very small. Nevertheless, we will neglect the convolution in Equation 
(10c), thus introducing a source of error in calculating R(R). It is usually 
the case that - F g = F^, where the subscript r refers to the lidar 
receiver. Often the gain of the receiver is adjusted between firing at the 
atmosphere and firing at the calibration target due to the much larger signal 
from the target. In addition, the gains and G g may depend on the 
received powers P, and P and on the bandwidth of the received power. 

D S 

Obviously, the bandwidth of P y (t) will be much greater than the bandwidth of 
P^(t). These effects on the receiver gain can introduce errors if not 
accounted for. We will Ignore these effects in what follows and simplify 
Equations (10b) and (10c) to 


V b (t) - G b , F r {P b (c» 

(12a) 

v s ( t ) = G s x F r (P s (t)} 

(12b) 


Combining Equations (4), (6), (11), and (12), we find that the first four 
terms in the right side of Equation (7) become 


P 

E 


tb 





0 


{v b (t)/ y 

T 


vP tb (t)dt 


0 



( t )dt 




F r 1{V s (t)/G s }dt 


2R 

s 


(13) 


c 


where the factors G^ in Equation (11) have canceled and {x} repre- 

sents the inverse operation of F^{x}. For example, the inverse of /x and 
2 x 

Anx would be x and e , respectively. At this point we observe a major 
difference between the linear, square-root, and logarithmic lidar receiver 
functions F . As stated earlier, the gains G, and G are often purposely 

IT D S 

made different due to the larger return from calibration targets than from the 

atmospheric aerosol. In practice, it is much easier to know the ratio of 

these gains ^ g /G^ than it is to know both their values independently. From 

Equation (13) it is clear that the ratio G /G^ is sufficient to calculate 

$(R) for linear and square-root lidar receivers, where F~*{x} s x and 

x , respectively. However, for a logarithmic receiver, F^ {x} = e' , 

and the actual values of G, and G are required to solve for P ( R ) - This 

b s 

suggests that logarithmic receivers are not appropriate if lidar aerosol 
backscatter data are to be calibrated using hard targets. 

It has been assumed in this section that a sufficiently fast 
transient digitizer is used to record C t ) » V fcs (t), V^(t), and V^(t). 

Another type of A/D converter that is often used is an integrating A/D con- 
verter. These devices usually have a variable gate width and integrate the 
voltage waveform during the gate period. It is clear from Equations (4) and 

(11) that an integrating A/D converter is acceptable to record the pulse 
energy signals, provided that the detector is linear in power. From Equations 

(12) and (13), however, it is apparent that the use of an integrating A/D 

converter would be appropriate for V^t) and V g (t) only if the function 

F^ is linear in power. For nonlinear lidar receivers, errors would be 

caused when an integrating A/D was used, especially with the higher-frequency 

calibration-target return signal V (t). 

s 
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SECTION 4 


DATA PROCESSING ERRORS 

4.1 SIGNAL-AVERAGING 

We have derived an expression for P(R), given by Equations (7) 
and (13), for the case of a single lidnr pulse directed into the atmosphere 
and a single lidar pulse directed at a calibration target. The energy of each 
pulse was used to normalize the return signals. It was assumed that the 
slgnal-to-noise ratios were adequate and that both speckle and atmospheric 
turbulence effects were not present. In reality, the received-power profiles 
from both the calibration target and the atmospheric aerosol will fluctuate 
from pulse to pulse due to speckle and atmospheric turbulence. These sources 
of signal fluctuation and the desire to improve low signal-to-noise ration 
usually make it mandatory to signal-average the results of many Ildar pulses. 

We now consider a measurement consisting of N^ pulses fired into 

the atmosphere and N ( pulses fired at the calibration target. We assume all 

s 

the terms in Equation (7) remain constant during the measurement except the 
four terms given in Equation (13). (Of course, the extinction profiles 
cij (R) and a^R), as well as B(R) itself, may change during the measurement 
time, due to temporal changes of the atmosphere or due to spatial movement of 
the lidar system, thus introducing a source of error. This again causes a bias 
against low B values that may exist for narrow spatial or temporal ranges.) 
Since normalization by the pulse energy will be discussed in a Later section, 
it is assumed here that all pulses are Identical, allowing cancellation of the 
two pulse-energy integrals in Equation (13). Under these assumptions, Equation 
(13) becomes 


N, 


OP’^’r'AL, PAGE JSf 
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<p b (t)> 

<1 > 
s 




’b 1-1 


2R 

— — ® -f T 
N c p 

I- E I 'rX^V dt 

S J- L 2R 


( 14 ) 


where < > indicates an ensemble-average and the superscripts i and j on v b (0 

and V^(t) represent the pulse-to-pulse variation in the recorded waveforms. 

Often, however, Equation (14) is not used to process signal-averaged lidar 

1 i 

data. Instead, the recorded voltage waveforms V b (t) an< ^ V-(t) arc? 
averaged and only afterward is the inverse operator applied. This 
often used but incorrect expression Ls 


-i| l 


<P b (t)> 
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b b 1=1 


vjct)! 
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s s 


N 
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V J g (t) 


dt 


(15) 


If F^ is a linear function, then Equation (15) is identical to Equation (14). 
It is interesting to calculate the error caused by using Equation (15) under 
typical lidar aerosol and target backscatter conditions, when F^ is a 
nonlinear function such as /P or £nP. 

Let X be an independent random variable with a probability 
density function p(x) and mean value <x>. Let Y = (1/N) [Xj + X2 + ... + X^] 
have probability density p(y). By the central limit theorem, for large N, 
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p(y) will be a Gaussian density function with the mean value <y> « <x>. This 
•naan value is given by 


<y> 


i 

/ 


<x> * J x p(x) dx 


(16) 


and represents the correct signal-averaged result given by Equation ( 1 A ) . We 
now consider a function of x, F^(x). If we signal-average F^(x), we again 
obtain a Gaussian density function with a mean value equal to <F [ .(x)>. 
Following the incorrect technique of Equation (15), we now perform the inverse 
operation: 


F r 1{<F r (x)>} 




(x) p(x) dx 


(17) 


The ratio of Equation (17) to Equation (16) yields the error caused by the 

signal-averaging technique of Equation (15), and this error obviously depends 

on p(x) and F (x). 

r 

We may approximate the probability density function of the returned 

power in a typical lidar system due to speckle and atmospheric turbulence by 

the Gamma density function of order M. This function is shown inset in Figure 

2 and ranges from an exponential for M = 1 to a Gaussian-like function (but 

with x > 0 only) for large M. The exponential (M = 1 ) form is predicted for 

fully developed speckle with no aperture-averaging. The ratio of Equation 

(17) to Equation (16) has been calculated for the Gamma density function and 

for F {x} == / x and Anx. For F {x} = / x the ratio is 
r r 


</ x > 2 
<x> 


TT 

M 


1 x 3 x 5 x 


x (2M - 


(M - 1)! x 2 


M 


1 ) 


(18) 
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Figure 2, Calculated error in the estimated mean value of a 
random variable x due to averaging the functions 
/ x and £nx vs. M where the probability density 
function of x is the Gamma density function of 
order M. 
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which approaches 1 for large M. The percentage error is plotted In Figure 2 
and reaches a maximum of -22% when M * 1. Referring to Equation (7) wo see 
that fl will be underestimated by 22?,, When F {x} « i’nx, the ratio in 
given by 


aR ... L<»^>) . . (0.561) M" 1 


r M-l . “i 

pxpE ~ 

n-1 n 


(19) 


which also approaches 1 for large M. As seen in Figure 2, the error in this 
case is greater than for F {x} “ /x and reaches a maximum of -44% for M * 1. 
Ernie (1975) addressed the closely related problem of comparing /X ("linear") 
and AnX (logarithmic) radar receivers with X ("quadratic") receivers. He 
treated only the M * 1 exponential power distribution (Rayleigh amplitude 
distribution) case and solved for both the error and standard deviation of the 
power estimate as a function of the number of samples, N. For large N his 
results are identical to our M = 1 results. He found that for large N, a 
"linear" receiver requires 1.09 times more samples than a "quadratic" receiver 
to achieve the same standard deviation, and the logarithmic receiver requires 
1.64 times more samples. 

It is clear that the signal-averaging technique of Equation (15) is 
unacceptable when the receiver function F is nonlinear, and -he correct 
technique of Equation (14) should be used. Only if p(x) and F^{x} are well- 
known could the expressions of Equations (16) and (17) be used to calculate a 
"correction" factor for the derived values of 8(R). 

4.2 PULSE-ENERGY NORMALIZATION 

In the last section it was assumed that all the laser pulses had 
Identical energies and, therefore, the terms E^ and E tg were equal and 
canceled in Equation (7). In actual practice, the transmitted pulse energy 
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will vary from pulse to pulse and proper normalisation by pulse energy is 
required to avoid errors. Expanding Equation (13) to Include correct signal- 
averaging, we obtain 



N b F r 1 {v b (t)/G b } 
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where the pulse indices 1 and j are now also necessary on the pulse, duration 
t for each shot, and the factors G fc in Equation (11) again cancel. If 
Equation (20) is inserted into Equation (7), we obtain the proper pulse-energy 
normalized, signal-averaged expression for <f?(R)>. 

Often, however, normalization by the pulse energy Is not done on a 
pulse-by-pulse basis, but rather the signal voltages and pulse energies are 
summed separately and the sums are ulvlded at the end of the measurement. It 
is interesting to calculate the resultant errors due to this technique. We 
may generalize either the numerator or denominator of Equation (20) for this 
analysis and use a more convenient notation. Let the desired ensemble-average, 
including the pulse-by-pulse normalization, be given by 
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UD 


where 1 b the return "signal" and E^ Is the pulse energy of the ith pulse. 
Let S' oe the incorrectly calculated ensemble-average as described above: 


S' 



1 N 
i-1 


( 22 ) 


If all the pulse energies E^ equal Eq, a constant, then S' - S, and either 
technique Is satisfactory. We now assume that E^ - Eq + e^, where Eq is the 
average pulse energy and represents the deviation of the ith pulse energy 
from Eq. Since it is obvious that S ' * S for large pulse-energy fluctuations, 
we will assume the deviations are small, l.e., << Eq. If we insert the 

expression for E^ into Equation (21) we obtain 
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ana inserting it into Equation (22) gives 
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Combining Equations (23) and (24) and using the fact that e^/ E q « 1 and 
that <e - 0, the fractional error is given approximately by 


Sj_-S 

S 


BB 



( 25 ) 


Obviously, the error goes to 0 as e^/Eq approaches 0. The si 2 e of the error 

will Increase to the extent that K. and e J are correlated. However, speckle 

i i 

an 1 atmospheric turbulence will tend to reduce this correlation. Nevertheless, 

this source of error should be avoided by normalizing to pulse energy on a 

pulse-by-pulse basis. 

4.3 RECTANGULAR PULSE ASSUMPTION 

In Section 2.1 it was shown that Equation (2) was an exact expression 

for the aerosol backscatter power profile, P^(t), provided that the aerosol 

backscatter coefficient, f<(R), and the telescope overlap function, 0(R), were 

slowly varying with respect to the spatial distance ct^/2. In going from 

-2 

Equation (2) to Equation (3), both the R and exponential terms were removed 
from the integrand in order that the remaining integral would reduce to the 
transmitted pulse energy, an easily measured quantity. All terms removed from the 
Integral in Equation (1) were assigned their values at the midpoint of the 
integration range, R^. The rationale for this simplification often starts with 
the assumption that the transmitted pulse, p fc ^(t), is rectangular with power Pq 
[W] , duration t , and therefore energy PqT^. Thus P fc ^ is removed 
immediately from the integral in Equation (2) and replaced with the constant Pq. 

At this point the range of integration is assumed to be narrow enough that the 
values at the midpoint, R^, are used for the remainder of the integrand. 
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A careful look at Equation (2), however, reveals that the R ' and 
exponential terms in the integrand combine to form a weighting function for the 
laser pulse profile P^Ct). This function favors the tail of the laser pulse 
(smaller R) in a complex way that depends on R, a^CR), and r . We have 
evaluated the percentage error in Equation (2) that results from assuming that the 
pulse shape is rectangular with power Pq when the actual pulse shape consists of 
two rectangular sections of unequal power and duration, as shown in Figure 3. 
Keeping the total pulse duration constant at t , we let the initial pulse power 
be Pj £W] for x^/a seconds and then provide for a "tail" with power P^/b 
for the remainder of the pulse. This bi-level pulse profile is a good 
approximation to the output profile of many pulsed lasers that exhibit an initial 
gain-switched spike followed by a lower-power, longer-lasting tail. The pulse 
energy is fixed at PqT by letting 


> \ 

1 ya + b - 1 J 


If a - b = 1 , Equation (26) reduces to the rectangular pulse. Figure 3 shows 

the resultant percentage error vs. range for several values of the parameters 

(a,b), for the total pulse duration t , and for a constant attenuation 

coefficient a^(R) = 0.125 km . It is clear that the percentage error can 

be very large, especially for large values of (a.b), that the error increases 

with increasing pulse duration, and that the error generally decreases with 

increasing range. Calculations using other values for show the 4- the 

percentage error is not strongly dependent on a^. For example, with t “ 5 ps, 

t = 4 y s, and (a,b) (8,8), where the error is large, letting a = 0.4, 

V b 

0.125, and 0 km yielded percentage errors of 72%, 68%, and 66%, respectively. 


23 


% ERROR 








These results show that F(R) will be overestimated by a range-dependent 
factor if Equation (3) is used. Only by using Equation (2), or correcting by 
the range-dependent factor, can this source of error be eliminated. Even if 
Equation (2) is used to solve for fl(R) by determining the pulse profile 
P tb (t ) and the atmospheric extinction profile a^(R), the return power 
profile ( t ) is not linearly proportional to pulse energy, making 
normalization to pulse energy very difficult. Only if the transmitted-pulse 
profile of each pulse were constant, so that P^Ct) could be written as the 
product of a pulse-energy term and a pulse-profile term, could Equation (2) be 
used to eliminate this "pulse-profile" source of error. Of course, both g(R) 
and 0 (R) may be varying significantly over distances of 0x^/2 and thus 
Equation (1) should be used, making the determination of P(R) very difficult, 
since the profiles f(R) and a^(R) are both unknowns in the measurement. 


SECTION 5 


TELESCOPE OVERLAP FUNCTION 

The range-dependent telescope overlap function, 0(R), was 
introduced in Equations (1) and (5) as an important term that affects the 
received-power profile from both the atmospheric aerosol and hard targets. 
Essentially, 0(R) is defined as the fraction of the transmitted pulse energy 
that is within the receiver's f ield-of-view. Obviously, factors that must be 
considered are the physical separation of the transmitting and receiving 
telescopes, their co-alignment or misalignment, the range-dependent spatial 
power profile of the transmitted beam, and the effective range-dependent 
spatial profile of the receiving telescope. Pulse-to-pulse variations in the 
transmitted spatial power profile would therefore produce pulse-to-pulse 
fluctuations in the telescope overlap function. The effective spatial profile 
of the receiver is due to many factors, including a central obstruction in the 
telescope, the finite size and shape of the detector element, sensitivity 
variations over the detector area, the variation in the position of the focal 
plane with range, and diffraction. For systems employing heterodyne 
detection, the local oscillator's spatial profile, polarization, and alignment 
at the detector plane must be considered. 

Often the telescope overlap function is neglected in the lidar 
equation, especially when coaxial lidar systems are being analyzed. However, 
all of the factors contributing to 0(R) for side-by-side lidar systems apply 
to coaxial systems except that the physical separation of transmitter and 
receiver is very small. Sassen and Dodd (1982) have analyzed the behavior of 
0(R) for both Gaussian and uniform transmitted-pulse profiles, and for various 
values of transmitter and receiver divergence, and transmitter-receiver 
misalignment. They have shown that 0(R) is strongly dependent on alignment, 


especially for typical Ildar systems, which employ narrow transmitter and 
receiver beamwidths. Referring to Equation (1), we see that experimentally 
determining 0(R) at a given wavelength requires that a very short rectangular 
laser pulse be fired into an atmosphere with well-known and slowly varying 
profiles 8(R) and a^(R). Since this experiment would be difficult, a 
model tor 0(R) is usually used instead. 

We have modeled the overlap function for the side-by-side telescope 
lidar system that is being used at JPL to measure the atmospheric aerosol 
backscatter coefficient in the troposphere in the 9 to 11 pm region (Kavaya 
at al., 1983). Figure 4 shows the modeled overlap function vs. range at the 
i(JP(.ZU) C0„ laser wavelength (10.39 pm). The indicated parameters (X,Y) 
represent the position (misalignment) of the detector element in thousandths 
of an inch (mils) with respect to the optical axis, which is parallel to the 
transmitted pulse, and lying in the focal plane (R = °°) of the receiver. 
Therefore, (X,Y) = (0,0) represents exact centering of the detector on the 
optical axis and (X,Y) = (-3,0) represents a 5-mil (127-pm) movement of the 
detector from the centered position. A 3-mil displacement is equivalent to a 
U.U74-mrad angular misalignment. The centers of the 15.5-cm diameter trans- 
mitter and receiver telescopes in our lidar system are displaced by 25.5 cm in 
the "X" direction and 7.3 cm in the "Y" direction. For our sign convention, 
negative values of X result in the crossing of the transmitter and receiver 
optical axes at some positive range from the telescopes. This range is 
indicated in Figure 4 for the appropriate curves. Note that the peak of these 
"negative X" curves does not occur at the crossover range but that the largest 
overlap function value at any range is given by the curve corresponding to 
that crossover range. The finite size (300 pm square) of the detector is 
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accounted for by averaging the center and four corner positions of the 
detector. In addition to assuming uniform detector responsiveness and using 
the five-point average, we make other assumptions in the model, including that 
there is a uniform local- oscillator , negligible receiver focal plane variation 
with the range beyond = 1.1 km, and a uniform, circular, transmitted- 

pulse spatial profile that diffracts at 150% of the diffraction limit. Al- 
though these assumptions limit the accuracy of the model, experimental results 
obtained by firing the Ildar system with various values of (X,Y) agree quite 
well with the curves in Figure 4 and also allow determination of the optimum 
(0,0) detector position. 

It is clear from Figure 4 that the overlap function varies signi- 
ficantly with very small detector misalignments. Not only is it important to 
know 0(R) when reducing the aerosol backscatter power profile, but it is 
especially important when using a hard target for calibration. A common range 
used with our calibration target is 2 km, which is indicated in Figure 4. The 
calculation of 3 at any range R will involve the ratio of the telescope over- 
lap function at the target's range, R , to its value at R as is shown in 
Equation (7). It is clear from the curves in Figure 4 that assuming an in- 
correct telescope overlap function (or neglecting it) can cause significant 
errors in calculating 3(R). For example, when the system is perfectly 
aligned, i.e., when (X,Y) = (0,0), then the ratio of the overlap function's 
value at 16 km to its value at 2 km is approximately 1. However, if the 
system were misaligned by only 220 prad, e.g., (X,Y) = (-15,0), then the 
same ratio is approximately 6. If correct alignment was mistakenly assumed, 
or if the telescope overlap function was neglected, then a ratio of 1 would be 
used and the value of 3(16 km) would be overestimated by a factor of 6. 


An additional problem in determining 0(R) occurs with airborne Ildar 
systems. This is due to the fact that once the laser pulse is transmitted, 
the receiver telescope will move with respect to its 1*0 position. Since 
calibration of the system with a hard target is most convenient on the ground 
before takeoff, the overlap function for the stationary lidar system will not 
apply to the moving (airborne) Ildar system. Fortunately, the sensitivity of 
0(R) to the displacement of the receiver will be greatest at short ranges, 
where the displacement will be the least. Whether the system is airborne or 
not, the telescope overlap function of every Ildar system should be modeled as 
accurately as possible at each wavelength of interest and checked 
experimentally. 
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SECTION 6 


HARD-TARGET CALIBRATION 

The hard-target parameter p* was Introduced In Equation (5) as 
the power reflected from the hard target per steradiun towards the receiver, 
divided by the Incident power. As is clear from Equation (7), a quantitative 
evaluation of this parameter is required to calculate P(R) when hard-target 
measurements are being used for calibration. It has been shown (Kavaya 
et al., 1983) that the evaluation of p* for a given calibration target, 
lidar wavelength, and lidar geometry is not trivial. Although hard targets 
are frequently used to calibrate lidar systems, there is a surprising lack of 
published information about their reflectance properties. The derivation of 
p * and a technique for its evaluation will be outlined in this section. 

The geometric reflectance properties of a flat, uniform, isotropic 
surface can be described by the bidirectional reflectance-distribution 
function (BRDF) [sr *3 defined by Nicodemus et al. (1977) as 
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In Equation (27) is the reflected radiance [W m sr ] in the direction 

_2 

(0 ,<(> ) due to the incident irradiance dE » L. cosQ. dm [W m ] confined 
r r iiii 

to the solid-angle element dm, in the direction (0. ,d>.); 0 is the polar angle 

x i i 

between the ray and the average surface normal; <t> is the azimuthal angle 
between the ray's projection in the surface and a reference direction in the 
surface; and the subscripts i and r on L, E, 0, and (j> refer to the incident 
and reflected radiation, respectively. The radiance In the direction (0,<f>) 
is defined by 


L(0,*) 


d 2 <t> 

dA cosQ dto 
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where cl <t> is the element of radiant flux [W] through the element of area 
dA in the direction (0,*) and within the solid-angle element dm. The 
element of the projected area perpendicular to the ray direction (0,<i>) is 
dA cos0. Equation (27) does not treat Interference, diffraction, transmis- 
sion, absorption, fluorescence, or polarization effects, and it is assumed 
that the illumination is monochromatic, uniform, and Isotropic. The BRDF is 
an unmeasurable derivative capable of values from 0 to infinity. Real 
measurements always involve an average of f^ over finite intervals, e.g., 

Am and AX. 

Nicodemus et al. (1977) define nine reflectances, allowing for 
directional, conical, or hemispherical geometry of the Illuminating and 
reflected radiation. The most basic quantity, the biconical reflectance, is 
defined as 


p(« i ;m r ) - 


f / f (0.j<j>.;0 ><f> ) cosO, cosf) dm, dm 
w r W 1 r 1 1 r r i r i i 

cos&, dm 


(29) 


where p(m^;m r ) is the ratio of the radiant flux in the direction (0 r »<f> r ) 

within m to the incident flux in the direction (0.,<J). ) within m . (We 
r ill 

will avoid the use of the projected solid angle dfi = cosod'o as it is 

nonphysical and can lead to confusion.) A Lambertian (diffuse) surface is one 

for which the reflected radiance is isotropic (L^ is constant, independent of 

0 and <j> ) regardless of how it is irradiated, and thus its BRDF, f , 
r it it j d 

is necessarily a constant. From Equation (29) we see that 
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If the solid angle to confining the illuminating radiation is 
small enough that we can consider to bo constant over the limits of 
integration, Equation (29) reduces to 


p(0 ,h> jot ) » / f (o ,rf> ;0 ,<j> ) cosn dw 

i 1 t J « r r i’ 1 r' r c r 


( 31 ) 


which is the directional-conica L reflectance, l.e., the illuminating radiation 

is from a specific direction, and the reflected radiation is measured in a 

cone of solid angle w . The reflectance measured by an ideal integrating 

sphere, p TC , is given by Equation (31) with w * 2n , 
is r 

2tr tt / 2 

P IS S p ^ 0 i , *i* 2ir ^ *00 f r^ < W 0 r ,<! V cosn r sin °r d °r dA r (32 ) 


which is the directional-hemispherical reflectance. If the small assump- 
tion is not valid, Equation (29) should be used for p TC with co “ 2ir. If 
the monochromatic illumination assumption is not true, an average over AX will 
be measured. The common term albedo is the hemispherical-hemispherical 
reflectance, l.e., p(2tt;2tt). 

In Ildar applications both the illuminating and receiving solid 

angles are very small. For example, a 15-cm diameter telescope subtends only 

—8 

7.1 x 10 sr at a range of 500 m. Thus we may further simplify Equation 
(31) by letting + 0, yielding 


P(0, ><!> . SO ><f> ) “ <»> <f (A. ><\>, »©. >0> cosR (33) 

II rr rriirr r 

where <f r > is the average value of f over the finite solid-angle intervals 

w . and w . 
i r 
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From Equation (33) It may be concluded that the received power from 

a Lambertian target (f r * constant) will vary as cosO , However, this 

conclusion is not always correct. It is important to consider the relative 

sizes of the target, the illuminated area on the target, and the area of the 

target viewed by the detector, os shown in Figure 5. 

Recall that in Equation (33), p is the ratio of the reflected 

flux to the Incident flux. The incident flux is not necessarily constant but 

is given by the product of the power per illuminated target area and the 

viewed illuminated target area. Figure 5a depicts the similar cases A^ < 

A <A or A, < A <A, where A, "A, /cos©, is the area illuminated 
rs i s r llo 1 

by a laser with constant power P, A^ e A^/cosO^ is the area viewed by the 

detector, and A, is the target area. Note that the illuminated spot is 
s 

smallest In this case, and, therefore, may not approach 90°. The power 

per illuminated target area is given by PcosO^/A.^, and the viewed illuminated 

target area Is given by A /cosfl^. Thus the cos©^ terms will cancel, leaving 

only the cos© dependence of p in the expression for received power. In 
r 

Figure 5b we have the cases A < A, < A or A < A < A, ; the receiver field- 

of-view at the target is smallest. Here the power per illuminated target area 

is again given by PcosQ^/A^, but the viewed illuminated target area is now 

given by A^/cosO . When these angular dependences are combined with Equation 

(33), we see that the overall dependence of the received power is given by 

cos©^. As seen in Figure 5, when a colllnear backscatter (e.g., lidar) geometry 

is considered with ©^ 55 ©^ = 0, all the above cases exhibit a cos© dependence. 

Finally, Figure 5c depicts the cases A < A < A or A < A < A , ; the target 

s i it s r x 

area is smallest. Once again the power per illuminated target area is 

Pcos0./A, , but the viewed illuminated target area is just A . When these 
l io s 
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Figure 5. Pictorial representation of three possible reflectance 
geometries: (a) the Illumination area at the target is 

smaller than both the target and the detector field-of- 
view at the target, (b) the fleld-of-view of the detector 
is smallest, and (c) the target area is smallest. 
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factors are combined with Equation (33), we obtain a cosO^cosQj, depend" 

ence. In this case the collinear geometry yields a (cosO) dependence. 

Of course, a more complicated geometry than those cases in Figure 5 is 

possible, e.g. , if , A r , and A g have boundaries that cross each other. 

It is clear that the geometry of any reflectance experiment is very important, 

especially when laboratory target reflectance data are used to characterize 

targets, which are in turn placed in a lidar geometry for calibrating lidar 

systems. Seldom in the literature is the geometry of a reflectance experiment 

defined in terms of the case? of Figure 5. The Lambertian reflectance behavior 

is predominantly given as cosO r (small A^), although examples of cosQ^cos©^ 

(small A ) can be found as well as examples of cos©, (small A ), especially in 
s i r 

astronomy. 

Later in this section we will discuss target measurements using a 

cw CO 2 laser in a collinear backscatter geometry (©^ = © r = ©), where 

A ~ A < A . These measurements provide collinear 0 dependence of each 
i r s * 

target's reflectance, showing various deviations from Lambertian (cos©) 
dependence. It is tempting to assume that each target's © dependence curve 
could be scaled absolutely by using integrating-sphere data and the area under 
the © dependence curve. However, this is not possible for surfaces with a 
general BRDF. The integrating-sphere data depend on the BRDF, as shown in 
Equation (32), while the backscatter data will depend on 0^, as shown in 
Equation (33): 
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where primes are used to differentiate from the angles in the integrating 

sphere case. We may let 0' =0' =0 and <b ! = <p ' = 0., where we use the 
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eollinear geometry and the assumed Isotropic property of the target (f 
independent of <|> ). If we now multiply S'"} by sinO and integrate over 
0, we obtain 

tt/2 tt / 2 

X * / S(0) sinO d0 * w' / <f (0,A ;0,<f> )> cosfi sinO do (35) 

0 r 0 r l l 

In comparing Equations (32) and (35) we see that we may not use Pj S to 
le S(0) unless certain conditions hold. First, the backscatter experi- 
ment's solid angles must be small enough to allow the assumption that <f r > = 
£^. Second, f^ mu. t. be independent of to allow the parameter A 
of in Equation (32) to be set equal to <i>^. Finally, f^ must be 
independent of 0^ to allow the parameter 0^ in Equation (32) to be set 
equal to 0^. Under these conditions the integrals of Equations (32) and (35) 
are proportional to each other. However, the requirement that f ^(O.^ 
be independent of 0^, d>^, and ^ is nearly as restrictive as the Lambertian 
assumption of constant f^ (i.e., independent of all four parameters). 

It is common to find expressions in the literature that give the 
parameter p* as either p, p/n, or (p cosO)/7r, where p is called the 
"reflectivity" of the target. If the BRDF of a calibration target is 
sufficiently known, we may more exactly define p*. Under the assumptions 
leading to Equation (29), e.g., uniform and isotropic illumination and target 
surface, we see that the general expression for p* in Equations (5) and (7) 
is 


/ w L. W WV C0Sfl i C0SPl r d( °i dW r 

p*(u) . ;u> ) = — - 

i r t 

to J cos 0 . dw . 
r J to, i \ 
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If small solid angles can be assumed, Equation (36) becomes 


P*^, 1 ^; 0 ..,^) “ <f r C° i »^ ;L > n r » f t> r )> cos^ r 


and if colllnear geometry Is used, Equation (37) becomes 


( 37 ) 


P*(o,<f>;0,'!>) = <f (0,*;0,r|>)> cosn (38) 

2 

If the geometry of Figure 5c applies, cosO is replaced by (cosO) in Equation 
(38). If the target is Lambertian (diffuse), then Equation (38) becomes 


P(w . 5 2tt ) p 

p*(0,fj);Q,(j)) = cos0 = — cosO (39) 

It is apparent that considerable calibration error may result from replacing 
p* in Equation (5) or (7) with either p, p/t, or even (p cosO)/-* when the 
experimental geometry resulting in the reported value of p is unknown; when 
the center wavelength, bandwidth, and polarization of the illumination may 
have differed from that in use; and when the target surface is most likely not 
Lambertian or reproducible. 

The ideal calibration target would have a well-known BRDF and would 
be inexpensive, easy to fabricate, durable, and reproducible. It has been 
shown that quantitative characterization of a target using integrating-sphere 
data and laboratory backscatter 0 dependence data is only possible if 
f (0 ,<f>.;0 ,<j) ) is independent of 0,, <*> . , and 0 , a very restrictive condition 

X X X IT IT IX l 

that is satisfied by Lambertian surfaces. Thus, a Lambertian target possessing 
the qualities listed above would greatly facilitate the calibration process. 
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If the reflected radiation of a Lambertian target were measured using a 
collinear backscatter geometry (0 * 0^= O) with A^ ~ < A g , then 

Equation (39) predicts a cos0 signal dependence as 0 is varied. The lack 
of a cosO signal dependence would indicate that a hard target was not 
Lambertian, while the presence of a cosO dependence would strongly indicate 
that the target was Lambertian. We have fabricated a collinear backscatter 
apparatus at JPL, which uses a cw C0 2 laser for illumination wavelengths in 
the 9 to 11 pm region, to investigate candidate calibration target surfaces 
(Kavaya et al., 1983). The laboratory collinear backscatter apparatus is 
shown in Figure 6. The passively stabilized, linearly polarized, cw C0 2 
laser is grating-tunable in the 9 to 11 pm region. The laser radiation was 
chopped at 100 Hz and was directed onto both the target surface and a power 
meter with a ZnSe beam-splitter. The vertical polarization of the laser was 
perpendicular to the plane of incidence at the target. A 30-cm focal length, 
5-cm diameter BaF 2 lens was used to image the target surface onto an 
LN 2 ~cooIed, HgCdTe infrared detector. A preamplifier and lock-in amplifier 
(not shown) were then used to obtain the signal magnitude. A ZnSe polarizer 
that had a greater than 500:1 extinction ratio and 10-mm diameter aperture was 
positioned next to the detector and was used to select backscattered radiation 
with polarization either parallel or perpendicular to that of the illuminating 
beam. (The different reflection coefficients of the beam-splitter for the two 
directions of polarization were measured, and the data were appropriately 
scaled.) The target was carefully positioned to align the center of illumi- 
nation with the axis of rotation. The beam size at the target was ~5 mm in 
diameter and was smaller than the target even at 0 = 80°. Since backscatter 
from the power meter contributed a large background signal, the power meter was 
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removed during measurements, and the unblocked beam was allowed to travel 
several meters before striking anything. Background readings were then taken 
with a mirror that was placed between the beam-splitter and the target, and 
that also directed that beam across the laboratory. The background readings 
were found to be negligible. The data were normalized by the laser power. 

'ihe finite size of the lens limited the resolution in 0 to 3.2°. However, it 
did Help reduce the signal fluctuations due to speckle, with an estimated 100 
speckle lobes being averaged. Speckle effects were further reduced by averag- 
ing several measurements within ±0.5° of the desired value of 0. Later, a 
spinning target was Incorporated to essentially eliminate the speckle effects. 

We chose four target surfaces for this work: sublimed flowers of 

sulfur (S), flame-sprayed aluminum (FSA) , 20-grit sandblasted aluminum (SBA) , 
and 400-grit silicon carbide sandpaper (SC). The sulfur target was fabricated 
by mixing sublimed sulfur (,J.T. Baker 4088) with acetone to form a slurry 
(~ 1.4 g sultur to 1 ml acetone). The slurry was packed into a rectangular 
well that was machined in aluminum to a 4.4-mm depth. A microscope slide 
lubricated with acetone was used to trowel a uniform surface finish that 
appeared smooth to the naked eye. The target was then allowed to dry over- 
night. The flame-sprayed aluminum target was made by sandblasting 6061 
aluminum plate with 60-grit sand, priming with Ni alumina bonding agent, and 
finally flame-spraying with pure aluminum to a thickness of 0.25 to 0.38 mm. 
The measurements (Kavaya et al. , 1983) showed that the FSA, SBA, and SC tar- 
gets all deviated significantly from the Lambertian cos0 dependence, whereas 
the £ lowers - of-sulf ur target obeyed the cos© dependence within our experi- 
mental uncertainty. The worst deviations occurred with the SC target, whose 
signal was fairly independent of 0, indicating that retroref lection was re- 
sponsible. 
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If we let the symbol s represent the backsc.attered radiation with 
polarization perpendicular to the illumination polarization direction, and p 
represent the backscattered component parallel to the illumination, then the 
depolarization ratio, defined as s/(p + s), varied from approximately 10% to 
30% for the four targets. The largest depolarization ratio corresponded to 
the flowers of sulfur, which is consistent with the idea that Lambertian 
behavior is caused by multiple subsurface scattering. The preservation of 
polarization direction with specular reflection was supported by the fact that 
the FSA and SBA targets exhibited specular reflection components near A = 

0° only in their parallel backscatter profiles and not in their perpendicular 
backscatter profiles, and by the fact that the depolarization ratio of the 
"retroref lecting" SC target was the lowest (~10%). 

Although our measurements indicated that the f lowers-of-sulfur 
target was Lambertian at X = 9.6 and 10.6 um, and presumably also at other 
CO 2 laser wavelengths, it is not an ideal calibration target. Our lidar 
beamwidths and calibration target range (2 km) require a large target surface 
(~1.8 m square) that is durable enough for field use. The sulfur targets, 
however, are very fragile and inconvenient to make in large sizes. We have 
found that the smooth surface appearance is required for the Lambertian 
behavior. Also, little is known about the effects of aging or exposure on the 
BRDF of the sulfur. 

To take advantage of the Lambertian behavior of the sulfur, we have 
adopted a primary and secondary standard-calibration technique. A flowers-of- 
sulfur target, which exhibits Lambertian behavior, is made the primary 
standard. An integrating sphere is then used to measure the directional- 
hemispherical reflectance, P]- s > of the target as a function of wavelength 


with an Incident Illumination angle of 45° (Kavuya et al., 1983). Since 
Equation (39) applies to the sulfur target, p* may be calculated for any 
angle 0. A large, durable FSA target is then fabricated and made the 
secondary standard. This target is used in the field with Q * 45° as the 
lidar calibration target. Assuming homogeneity, a small piece of the FSA 
target, suitable for use in the apparatus of Figure 6, is cut off. 

Using the apparatus shown in Figure 6, a ratio of backscattar 
reflectance is determined between the primary sulfur target and the secondary 
FSA target. This ratio is measured at 0 - 45°, with the incident and 
detected directions of polarization matched to that of the lidar geometry, and 
at each wavelength of interest. Finally, the value of p* for sulfur, 
calculated from Equation (39), is multiplied by the measured reflectance ratio 
of FSA to sulfur, to obtain the p* for FSA that is Inserted into Equation 
(7). For example, we have used this technique to determine a value of 0.097 
for p* of our FSA lidar calibration target at the 10P(20) C0 ? laser wave- 
length. If FSA is assumed to be Lambertian, however, with a "reflectivity" of 
80%, which has been done in the literature, then the calculated value of p* 
would be 0.8 cos 45°/ 7r = 0.18. This value is 86% larger than our value. 

Ideally, the reflectance transfer from the primary to the secondary 
standard should be done under conditions that match both the integrating- 
sphere measurement and the lidar measurement as closely as possible. In 
practice, however, there are many departures from this ideal. The transfer 
angle, 0, should match the incident angle in the integrating-sphere 
measurements and the lidar-target polar angle. The solid angles m. and 
w in all three measurements should be as close as possible. (This may be 
impractical since lidar solid angles are so small. Thus the BRDF of the 


primary and secondary targets should not vary strongly with 0. This 
consideration makes retroref lecting targets a poor choice.) The sensitivity 
of the reflectance transfer ratio to and could be studied by varying 
the solid angles in the backscatter apparatus. A tunable cw CO 2 laser should 
serve as the illumination source for the integrating-sphere data with polariza- 
tion geometry identical to the other two measurements. (If the integrating- 
sphere measurements of sulfur are made with unpolarized light* then four 
measurements of t ic sulfur at 0 « 45° should be made with the collinear 
backscatter apparatus. A measurement should be taken with all four permuta- 
tions of parallel and perpendicular, incident and reflected radiation: ss, 

sp, pp, and ps. One-half of the sum of these four readings should then be 
ratioed with the single (e.g., ss) reading of the FSA target.) The same 
Lambertian sulfur target should be used for both the integrating-sphere and 
transfer measurements. The secondary target should be as homogeneous as 
possible so that the small piece used in the transfer measurement accurately 
represents the larger lidar target. (We have seen variations in the BRDF of 
different FSA targets that probably are due to the operator, the fabrication 
technique, the size of the target, etc.) In addition, the small piece cut 
from the secondary target should be kept with the larger target, and thus be 
subjected to the same environment, so that later transfer measurements can be 
used to update .the calibration and to study aging and exposure effects. 

It has been shown that the calibration-target parameter p*, which 
must be well-known for accurate calibration of aerosol backscatter data, de- 
pends intimately on the BRDF of the target surface, a function that can only 
De partially measured with a costly, time-consuming, and complex goniometric 
experimental measurement and that likely varies with the wavelength, polariza- 
tion, and bandwidth of the illuminating radiation, and with the subtended solid 


angles of the illumination and detector. Little is known about the variations 
of the BRDF of hard targets with these parameters, or the variations with time 
and exposure for a single target, or the variations among different samples of 
the same type of target surface. It appears at present that the best method 
for determining the value of p* for the lidar equation is to use the primary 
and secondary standard technique outlined above, Including consistent and coor 
dinated integrating-sphere, collinear backscatter reflectance ratio, and lidar 


measurements. 


SECTION 7 


ATMOSPHERIC MODELING 

The remaining two parameters of Equation (7) that must be known In 

order to calculate B(R) are a (R) and a (R), the atmospheric extinction 

s b 

profiles for the optical paths from the Ildar to the hard target and aerosol 

particles, respectively. Ideally, both profiles would be accurately measured 

at the time of each Ildar measurement and at the transmitted wavelength. 

Unfortunately, this Is very difficult to do. It Is clear from Equation (2) 

that the Ildar signal does not readily yield the atmospheric extinction 

profile when B(R) and n b ( R ) are both unknowns. An analytic relationship 

between B and is not possible at 10 pm since the extinction is 

mainly molecular while the backscatter is due to the aerosol particles. 

However, if the pulse duration is sufficiently short, if the overlap 

function Q(R) is well-known, and if both B(R) and a^(R) can be assumed 

to be independent of R, then the lidar signal in Equation (2) can be used to 

determine the extinction coefficient. These conditions are reasonable for 

horizontal paths if both 3 and n depend only on altitude. Thus the 

extinction profile ex (R) may be considered to be constant if a horizontal 

path between the lidar and the calibration target is employed. Furthermore, 

by aiming the lidar so that it just misses the target, and therefore obtaining 

backscatter from a much longer yet similar horizontal path, the value of 

a (R) = a may be determined. For a ground-based lidar, such as the 
s s 

coherent COg lidar at JPL, the horizontal path lies in the boundary-layer 

and a is called the boundary-layer total attenuation . 
s BL 

These steps are shown in Figure 7 in relation to the overall 
calibration process used at JPL to reduce the lidar backscatter data and 
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Figure 7. Flow diagram depicting the main components in the 
lidar calibration and 0 profile computation. 
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calculate vertical and horizontal fl profiles from vertical and horizontal 
aerosol backscatter data and from horizontal hard-target backscatter data. 

Note that the steps of data acquisition and data processing, as discussed in 
Sections 3 and 4, ore not portrayed in Figure 7, However, both the telescope 
overlap function and hard-target calibration, discussed in Sections 5 and 6, 
are shown in the calibration flow diagram. 

Obtaining the extinction profile ot ^ C R ) for the vertical aerosol 
backscatter data is much more difficult than obtaining a (R). Equation 

h 

(2) may not be used, since both a^CR) and fl(R) will depend on R. Direct 
measurement of extinction vs. height at the transmitted wavelength and at the 
time of each measurement would be very difficult and expensive. At CO^ 
laser wavelengths the total extinction consists of molecular extinction and 
aerosol extinction. If the CO 2 line is chosen to avoid resonances of trace 
species such as ozone, then the molecular extinction will be due primarily to 
CC >2 and water vapor. The extinction due to these two molecules may be 
calculated from altitude profiles of atmospheric temperature, pressure, and 
relative humidity. These profiles may be measured directly using rawinsonde 
or radiosonde ascent probes. An alternative to the Inconvenience and cost of 
ascent probes is to use a model for the atmospheric temperature, pressure, and 
relative humidity that is appropriate to the measurement location and season. 

As shown in Figure 7, this latter technique is used at JPL. The molecular 
attenuation vs. altitude above the boundary layer is calculated using the 
Midlatitude-Summer model profiles for temperature and water vapor (McClatchey 
et al., 1972), the CO 2 and H 2 O absorption-line parameters from the 1982 
AFGL atmospheric absorption-line parameter compilation (Rothman et al., 1983), 
and the water-vapor continuum absorption parameters given in LOWTRAN 5 
(Knelzys et al., 1980). 


48 


ORIGINAL PAGE 19 
OF POOR QUALITY 


* 


Figure 8 shows the calculated molecular attenuation for various 
CO,, laser lines* Within the boundary layer the total attenuation calculated 
t rom the horizontal aerosol backscatter is used, (A slant-path aerosol 
backscatter measurement through the boundary layer may also be employed to 
determine the height of the boundary layer.) The molecular (C0 2 + HgO) 
contribution to the boundary-layer attenuation is calculated from direct 
measurements of temperature, relative humidity, and pressure. This is 
subtracted from the total boundary-layer attenuation to obtain the aerosol 
contribution in the boundary layer. The aerosol attenuation normally 
decreases rapidly with increasing altitude above the boundary layer. A model 
for the aerosol altitude scaling factors (McClatchey and D'Agati, 1978) Is 
used to derive the aerosol contribution to attenuation vs. altitude. These 
altitude scaling factors are also shown in Figure 8. Finally, the total 
attenuation vs. altitude n^(R) (see Figure 7) is found by adding the 
aerosol and molecular contributions. 

Any errors in the two profiles, a (R) and a. (R), that are 

s b 

used in Equation (7) will obviously produce errors in the calculated profile 
3(R). If we let a s ( R ) and a ^( R ) represent the true atmospheric 
profiles, and a^(R) and a^(R) represent incorrect profiler that are 
used in Equation (7) to calculate B'(R), then the ratio of incorrect to 
correct 8 is given by 


B ' ( V 

exp ‘ 

R 

r s 

|_ 2 J «'(R')dR'| 

« R b> 

exp \ 

( R \ 

{ f S a (R')dR'> 



0 



« b (R')dR'J. 

a^(R')dR'}> 


0 


(40) 
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Figure 8. Molecular attenuation and aerosol extinction scaling factor 
profiles that are used above the boundary layer for com- 
putation of the total attenuation. A Midlatitude-Summer 
model of atmospheric pressure, temperature, and relative 
humidity .was used. The attenuation at 9.11 pm, which is 
not significantly affected by atmospheric carbon dioxide, 
is shown for comparison. 
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The ratio given by Equation (40) may be applied to the typical atmospheric 
modeling geometry shown in Figure l J, We let the lidar system fire into the 
atmosphere at zenith angle 0 from altitude Z^. The atmosphere is modeled 
by layers i that have total attenuation a^ and an upper altitude of Z^. 

When appropriate, the first layer is the boundary layer with attenuation a^ 

* and altitude = Z^. The horizontal path to the calibration 
target Is in the first layer, and the altitude Z at which 6 is being calcu- 
lated is assumed to lie in layer j, With this geometry, Equation (40) becomes 


S'( V 

m b > 


- exp {-2 Ma B ' L - « BL >} exp{- ^ [« BL V V'BL + WjT V]} 

* exp { - dW [“2 (Z 2‘ hO - “2 (Z 2 - x • • • 

X exp i - JL_. fa (Z. Z, ,.) - a: ,(z! Z' )1 ) 
l cos© [J" 1 J" 1 1“ 1 J" 1 1“ 2 J / 

x exp | - fa' Z ' - a. Z. + Z(a.- a()1) (41) 

l cos0 l J j-1 J j- * J 3 JJ- 


The effects of making specific errors in the atmospheric attenuation model may 
be easily derWed from Equation (41). For example, if the only error in the 
atmospheric model is in the boundary-layer total attenuation a^, then the 
ratio of Equation (41) for Z > Z^- L reduces to 


6, cy 

B <v 


exp 


{■ 


2(a' - 
BE 


a ) 
BL ' 


[R - 
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(Z BL - Z L )/cos0] 


(42) 


which is independent of Z (i.e., R^ ) . Note that a prudent choice of R s and 
0 in Equation (42) could eliminate errors in 3 due to any inaccuracy in 
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Figure 9. Pictorial representation of the layered atmospheric attenuation 
model. 
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a RL . As a second example, consider Che only model error to be the height of 
the boundary layer Z Then Equation (41) becomes 


r(V 


cosO (Z BL' Z BL )( “2“ “bI^} 


which is again independent of Z. If the bracketed exponents in Equations (42) 
and (43) are much smaller than 1, then the percentage error in P will be 
proportional to the exponents. 

We have investigated the altitude dependence of the error in 8 
that results from using various incorrect atmospheric attenuation models. 

Data from an actual backscatter profile taken on July 21, 1983, at 10.6 pm 
with D=0° were used in the investigation, and the results are shown in 
Figure 10. Our lidar altitude was approximately 0.4 km and the assumed correct 
atmospheric model was a Midlatitude-Summer model with a 1.5-km boundary-layer 
altitude and a 0.39 km * boundary-layer attenuation. Figure 10 gives the 
percentage error in P(Z) resulting from using four "incorrect" atmospheric 
models: 1) a Midlatitude-Winter model above the boundary layer, 2) a Tropical 

model above the boundary layer, 3) a Midlatitude-Summer model but with a 
boundary layer altitude of 1 km, and 4) a Midlatitude-Summer model with a 2-km 
boundary layer altitude. As predicted by Equation (43), the effects of using 
an incorrect value of are independent of altitude. The error due to 
assuming Zp^ = 2 km is ten times the error due to assuming Z ^ ~ 1 km. 

This is expected from Equation (43), since (o^ - a^) is much larger 
for the former case. The use of incorrect atmospheric models above the 
boundary layer results in an altitude-dependent error in B, as expected from 
Equation (41). In our example, the error increases with increasing altitude, 
reaching as much as 30% at 8 km. 
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Figure 10. Plot of error vs. altitude in calculated 3 values 

resulting from the assumption of incorrect atmospheric 
models or of boundary-layer altitude. The nominal back- 
scatter profile used actual July 21, 1983, backscatter data 
at the 10P(20) CO 2 laser line; a Midlatitude-Summer 
temperature, pressure, and humidity profile; and a 
boundary-layer altitude of 1.5 km. 
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It Is clear from these results that an incorrect atmospheric attenu 

ation model can easily lead to ±30% errors in the calculated values of 8. 

Correct modeling of the boundary layer is Important for upward-pointing lidars 

that pass through a significant boundary layer. Of course, the boundary-layer 

effect on downward-pointing Ildar data is much less significant. 

Other factors that may have to be considered in modeling the 

atmosphere, depending on the lidar geometry and operating wavelength, include 

multiple scattering (Tam, 1983), Rayleigh scattering (Sroga et al., 1983) and 

the effect of the laser pulse on the atmospheric parameters (Fowler, 1983). 

Multiple scattering and Rayleigh scattering become relatively more Important 

at shorter wavelengths and should be considered when using lidar at, e.g., 

2 

1 pm. Fowler (1983) indicates that transmitted pulse fluences of 1 J/cm 

or greater at CO 2 laser wavelengths may alter the characteristics of fog or 

aerosols composed largely of water. The pulse fluences of typical coherent 

2 

lidars used in atmospheric studies are below 0.1 J/cm . 
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SECTION 8 


CONCLUSIONS 

Several possible sources of systematic error in interpreting and 
calibrating lidar aerosol backscatter data have been reviewed in this report. 

In addition to the difficulties posed by the physics of the measurement, other 
areas of concern that were discussed include data acquisition, data processing, 
hard-target calibration of the data, and modeling of the telescope overlap 
function and the atmospheric attenuation. In each case, examples of the re- 
sulting error due to incorrect data reduction, data interpretation, calibra- 
tion, or modeling were given. The largest potential source of error was due 
to incorrectly modeling (or neglecting) the telescope overlap function. Al- 
though the other potential error sources were smaller, they each must be con- 
sidered in relation to the desired measurement accuracy of the backscatter 
coefficient . 

It is hoped that this discussion will contribute to more accurate 
and consistent calibration of lidar data through increased understanding of 
the necessary steps in data interpretation and calibration. Although the 
viewpoint pertained primarily to stationary, pulsed-CC^ lidar measurements, 
much of the analysis should be useful for ground-mobile and airborne lidars, 
for lidars operating at other wavelengths, and for cw-focused lidars. 

It is believed that careful consideration of the calibration- 
process steps described here will allow significant improvement in the 
accuracy, calibration, consistency, and inter-calibration between different 
systems of lidar aerosol backscatter data. 
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APPENDIX A 


SYMBOLS AND DEFINITIONS 


Symbol 

Definition 

Units 

A 

Area 

2 

m 

A 1 

Illuminated area 

2 

m 

A 

r 

Area viewed by receiver 

m 2 

A 

s 

Hard-target surface area 

m 2 

c 

Speed of light (2.998 x 10 8 ) 

m s" 1 

C 

Capacitance 

F 

E 

Irradiance 

„ " 2 
W m 

B 

Energy 

J 

E 0 

Average pulse energy 

J 

E i 

Energy of pulse i = Eq + 

J 

E tb 

Energy transmitted into atmosphere 

J 

E ts 

Energy transmitted toward hard target 

J 

f 

Frequency 

Hz , s' 

f 

r 

Bidirectional reflectance-distribution function 

sr" 1 

F 

Operator 

- 

F b 

Aerosol backscatter operator 

- 

F 

r 

Lidar receiver operator 

- 

F 

s 

Target backscatter operator 

- 

F t 

Transmitted-pulse operator 

- 

G 

Gain parameter 

- 

S 

Aerosol backscatter gain parameter 

- 
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Target backscatter gain parameter 
Transmitted pulse gain parameter 
Heaviside unit step function (0:t<0; l;t>0) 
Counting index 

Time integral of target return power 
Counting index 
Arbitrary signal 
Radiance 

Gamma density function parameter 
Number of events 

Number of aerosol backscatter pulses 

Number of hard-target backscatter pulses 

Telescope overlap function 

Polarization parallel to plane of incidence 

Probability density function 

Probability density function 

Power 

Power of rectangular pulse 
Aerosol backscatter received power 
Target backscatter received power 
Transmitted-pulse power 
Transmitted-pulse power into atmosphere 
Transmitted-pulse power toward hard target 
Resistance 
Range 

c[t/2 - t /4] 

P 


J 


W m~^ sr 


W 

W 

W 

W 

w 

w 

w 

9 

m 

m 


60 


m 


rt 

s 


t 


V 


V 


tb 


V 


ts 


X 

Y 


Z 


Z 


BL 


Z 


L 


u 


u b (R > 

a 

BL 

u s (R) 

6 

6(x) 


s 


1 


n 

o 


Range to hard target 

Polarization perpendicular to plane of incidence 

Arbitrary signal 

Time 

Voltage 

Aerosol backscatter recorded voltage 

Target backscatter recorded voltage 

Transmitted pulse recorded voltage 

Recorded voltage for transmitted pulse into 
atmosphere 

Recorded voltage for transmitted pulse toward 
hard target 

Random variable 

N _1 2X 

Altitude 

Altitude of top of boundary layer 
Altitude of Ildar apparatus 
Extinction coefficient 

Extinction coefficient along atmospheric path 
Extinction coefficient of boundary layer 
Extinction coefficient along target path 
Aerosol volume backscatter coefficient 
Impulse symbol 

Deviation of ith laser pulse energy from Eq 
O ptical efficiency 
Polar angle 


s 

V 

V 

V 

V 

V 

V 


x 

x 

m 


m 


tn 


-1 


m 


-I 


-1 


m 


-1 


-1 -1 
m sr 


x 

J 


-1 


rad 
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X 


ft 

p(w i Jw r ) 

P IS 

P* 


T 

T 

T 

T 

T 


b 

p 

s 

t 


4> 




0 ) 


Optical wavelength 
3. 14159 

Biconical reflectance 

Directional-hemispherical reflectance, pCO^,*^; 2n ) 
Target reflectance parameter 
Time constant 

Aerosol backscatter time constant 
Transmitted pulse duration 
Target backscatter time constant 
Transmitted pulse time constant 
Azimuthal angle 
Radiant flux 
Solid angle 

Incident irradlance solid angle 
Receiver solid angle 


m 


sr 

s 

s 

s 

s 

s 

rad 

W 

sr 

sr 

sr 


<x> Average value 
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